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ABSTRACT 

We present simulations of a model for the formation of a prominence 
condensation in a coronal loop. The key idea behind the model is that the 
spatial localization of loop heating near the chromosphere leads to a catastrophic 
cooling in the corona ( [Antiochos fc Klimchuk 1991|) . Using a new adaptive grid 
code, we simulate the complete growth of a condensation, and hnd that after 
~ 5, 000 s it reaches a quasi-steady state. We show that the size and the growth 
time of the condensation are in good agreement with data, and discuss the 
implications of the model for coronal heating and SOHO/TRACE observations. 

Subject headings: 
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1. Introduction 


Solar prominences and filaments are beautiful phenomena both in their appearance 
and in their physics. They are observed in Ho as large masses lO^^g of dense > 10^^ 
cm“^ chromospheric plasma which seem to float high up in the tenuous (~ 10® cm“^) 
solar atmosphere, the million-degree corona. They can have a myriad of shapes, but often 
resemble a long arched wall with widths ranging from 10® to 10® cm, heights ranging 
from 10^® cm down to the chromosphere itself, and lengths that can be of order the Sun’s 
circumference for high-latitude polar crown hlaments. Typical lifetimes of prominences can 
range from as short as hours for active region hlaments to as long as months for quiescent 
prominences. Prominences usually disappear by eruption. Large coronal mass ejections and 
two-ribbon hares are almost always accompanied by prominence eruption. Consequently, 
the structure of prominences, especially their magnetic structure, is vital for understanding 
the major solar drivers of geomagnetic disturbances. 


There are three central issues concerning the physics of prominences: How does mass 
condense in the corona, how is it supported, and why does it erupt? The coronal magnetic 
held undoubtably plays a vital role in all three issues {e.g. [Tandberg-Hanssen 1995|, Priestj 
19891). We believe, however, that these questions are essentially decoupled in that the 


eruption mechanism is due to the global topology of the held (Antiochos 1998, Antiochos^ 


DeVore, fc Klimchuk 1999| ), the support is due to the local geometry of a sheared 3D arcade 
( lAntiochos, Dahlburg, fc Klimchuk 1994| ), while the condensation formation is due to the 
properties of the coronal heating process and the cooling of coronal plasma by radiation 
and conduction ([Antiochos fc Klimchuk 1991| , Pahlburg, Antiochos, fc Klimchuk 1998|) . 

In this paper we focus on the third question of condensation formation and, in particular, 
address the dynamics of this process. 


When viewed with high spatial resolution, it can be seen that the mass of a prominence 
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is not one continuous structure, but is actually a collection of small condensations or knots 
that form on time scales of tens of minutes and have size scales of order a thousand km. 
( [Engvold 19761 , [Zirker 1989| ). Given typical coronal gas pressures, ~ 10“^ ergs/cm^ and 
magnetic held strengths 10 G, the plasma is low-beta and we expect that each condensation 
is isolated from the others by the magnetic held. Therefore, the canonical picture of a 
prominence is that of an arcade of hux tubes, or coronal loops, each of which has somewhere 
along its length a plug of cool material that is observed as an Ha knot, and the ensemble of 
all the knots forms the global structure that we call a prominence. 


An important point is that prominence knots cannot be considered as simply the 
condensation of coronal plasma because they generally contain more mass than the coronal 
portion of the hux tubes they occupy. Material must be brought up from the chromosphere 
in order to produce the large mass of a prominence, which has led many authors to consider 
the so-called siphon how model for condensation formation {e.g. pikehner 1971| , pZngvold"^ 


Jensen 1977|, Priest &: Smith 1979| , Pibes fc Unno 19^ Poland fc Mariska 1986|, |An, Bao, 


|fc Wu 19^ , |Mok et al 19901) . It has proved to be difficult, however, for such models to 
produce condensations of sufficient mass (Poland fc Mariska 1986|) . The observation that 
the knot contains a large mass implies a large amount of chromospheric evaporation {e.g. 
[Antiochos fc Sturrock 1978| ) and, hence, a large coronal heating rate in the loop. But a high 
heating rate implies high coronal-like temperatures, rather than the cool chromospheric 
temperatures needed for a prominence. 


We have described elsewhere how one can reconcile the seemingly contradictory 
requirements of high coronal heating and low coronal temperatures ( [Antiochos fc Klimchuk 
19911 , Pahlburg, Antiochos, &: Klimchuk 1998|) . The key idea is that if the heating in a 


coronal loop has a strong spatial dependence, then it is possible to have a high heating 
rate and still develop a massive condensation near the loop center. Note that if the heating 
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is uniform per unit volume or per unit mass, which are the usual assumptions in coronal 
loop models {e.g. [Rosner, Tucker, fc Vaiana 1978|, [Vesecky, Antiochos, fc Underwood| 
|1979|) , then the loop can only have a hot structure with coronal temperatures throughout. 
This is true even if the loop has a dipped geometry suitable for supporting prominence 
material ([Antiochos &: Klimchuk 1991|) . If the heating is sufficiently localized near the loop 
base, however, the plasma structure can change to one in which there are chromospheric 
temperatures in the dipped portion of the loop. 


In our previous paper ( [Antiochos fc Klimchuk 1991|) , we verihed this hypothesis with 
numerical simulation. We began the simulation with a coronal loop that had a dipped 
geometry and a spatially uniform heating. This loop developed a typical hot coronal 
equilibrium solution with a temperature maximum > 10®K at the loop midpoint and a thin 
transition region at the base. We then increased the heating, but only in a small region near 
the base. At hrst, this heating increase caused the temperature to rise throughout the loop, 
so that the central dipped portion became even hotter. Chromospheric evaporation was 
driven by the temperature rise, which caused the density and, consequently, the radiative 
losses to increase throughout the loop, even in the upper section where the heating was 
unchanged. Eventually, the increase in radiative losses produced a catastrophic cooling at 
the loop center, and the formation of a chromospheric condensation. 


It is worth noting that this temperature collapse is physically due to a loss-of- 
equilibrium rather than a thermal instability. For small increases of the heating near the 
base, it is possible to hnd static equilibrium solutions which have a hot corona throughout 
the loop. These solutions are characterised by a slight increase in the temperature at the 
loop midpoint as the heating increases. But if the base heating is increased sufficiently, 
then no static solution with a hot midpoint is possible, because that would imply too high 
a radiative loss rate there. The only static solutions are those in which the temperature in 
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the middle section of the loop is very low, T « 10^ K, so that the radiative losses there 
are reduced (Pahlburg, Antiochos, fc Klimchuk 1998|) . 

Although our previous simulation verihed the basic idea of the model, there was 
a problem with that calculation which prevented us from comparing the model with 
observations, especially the time scale for prominence formation. The code could not 
follow the condensation signihcantly into its growth phase. As the condensation formed, 
a transition region with very large temperature gradients formed on each side of the 
condensation, exactly like the transition region at the footpoints. This transition region 
moved down the loop as the condensation grew, but since our code used a hxed grid, it was 
not possible to resolve a moving transition region and, consequently, the simulation broke 
down almost immediately after the condensation hrst appeared ([Antiochos fc Klimchuk] 
|1991|) . In this paper we present simulations using a new code that has a fully adaptive grid 
( [MacNeice et al. 1998]) and, therefore, can follow a condensation throughout its evolution, 
allowing us to compare the growth times with data. These new calculations also include 
more realistic forms for the loop geometry and the heating. 


2. The Numerical Model 

Since the magnetic held in a prominence dominates the plasma, we make the usual 
assumption of a ID coronal loop model in which the held determines the loop geometry 
and all plasma quantities vary only along the held. Our code solves the standard set of 
transport equations for a ID plasma: 


Wt 


d d 

(pV) + ^^(pV^ + p) = -pat(s), 


( 1 ) 


( 2 ) 
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and 


+ ^{UV - 10-6t5/2^T) = -P^V + E{t, s) - n^k{T). 
at os os os 


( 3 ) 


In Eqnations (1) - (3), s is distance along the loop, i.e. along the magnetic flnx tnbe; p 
is the plasma mass density; V is the velocity along the loop; P is gas pressnre; g\\ is the 
component of gravity along the held; U = 3P/2 is the internal energy; is the 

Spitzer (1962) coefficient of thermal condnction; E{t,s) is the assnmed form for the coronal 
heating rate; n is the electron nnmber density; and A(T) is the radiative loss coefficient 
for optically thin emission {e.g. Cook et al 1989 ). Note that we assnme a fnlly ionized 
hydrogen plasma, so that p = 1.67 x 10“^'^ n and the eqnation of state is simply P = 2nkT. 


The code that we have developed ( [MacNeice fc Spicer 1998|) to solve these eqnations 
has two major improvements over onr previons code. First, in order to advance the solntion, 
it nses a second order Godnnov scheme with a MUSCL limiter ([Van Leer 1979|) applied to 
the characteristic variables. It is widely believed that higher order Godnnov schemes, which 
solve the nonlinear Riemann problem at each cell face, are the most robnst methods for 
nnmerically solving ID flnid flow problems like the one we expect for coronal loop plasma 
— a non-tnrbnlent compressible flow in which steep gradients and/or discontinnities are 
likely to develop ([Hirsch 1988] , [Golella 1990|) . 


Second, the code employs a fnlly adaptive grid by nsing the PARAMESH Adaptive 
Mesh Refinement package ( MacNeice et al. 1998|) . PARAMESH bnilds a hierarchy of 
snb-grids to cover the compntational domain, with spatial resolntion varying to satisfy the 
demands of the application. These snb-grid blocks form the nodes of a tree data-strnctnre. 
Each grid block has a logically cartesian mesh, and the index ranges are the same for every 
block. In this case we nse ID snb-grid blocks with 20 grid cells and 2 gnard cells at each end 
of the snb-grid. When a snb-grid block needs to be refined it spawns 2 child blocks. These 
two children have their own 20 cell snb-grids, and together cover the same coordinate range 















as their parent, but now with mesh spacing one-half that of the parent’s mesh spacing. 
Any or all of these children can themselves be rehned in the same manner. PARAMESH 
provides routines which manage the necessary communication between sub-grid blocks. 

An adaptive grid is essential for our simulation because we expect a thin transition 
region to form somewhere in the loop and to move with time. Unless this region is resolved, 
the plasma energetics and dynamics cannot be calculated accurately. Often, a tricky issue 
is the selection of a criterion for rehning and derehning the numerical grid. We find that 
density is the best quantity to use in determing the rehnement because it exhibits a strong 
spatial variation, both in the transition region where pressure is approximately constant, 
and in the chromosphere, where temperature is approximately constant. In the simulations 
presented here we use six levels of rehnement, which corresponds to a minimum grid spacing 
of ~ 5 km. This is sufficient to do an adequate job of resolving the transition region 
throughout the evolution. 

For these calculations we established a very simple rehnement/derehnement criterion. 
At each grid cell i we computed the error measure 

Ci = \pi - pi-i\/mm{pi,pi_i), (4) 

and for each sub-grid found the maximum value 

= max(eo,ei,... , 621 ). (5) 

Notice we include one guard cell at each end of the sub-grid (ie i = 0,i = 21) to enable 
rehnement in advance of an arriving feature of the solution. If emax exceeds 0.25, meaning 
that somewhere in the sub-grid block the density varies by at least 25% between grid cells, 
then that block is rehned. If emax is less than 0.05 on a pair of sibling sub-grid blocks, then 
the block is derehned provided the lower resolution representation of their parents does not 
exceed the rehnement threshold. 



The refinement process was permitted to range freely between refinement levels 6 and 
12. In this context, level 1 indicates that the entire computational domain is covered by 


1 sub-grid block. A level 1 sub-grid would have a cell size of Ax = L/20 where L is the 


total length of the computational domain. A level n sub-grid therefore has a cell size of 
2 -n+i^/ 2 o, The initial grid was uniformly refined at level 6, ie Ax = 2“®L/20 = 0.0015625L 


The sub-grids refined to level 12 have a cell size which is a factor 2® = 64 finer. 


The geometry that we assume for the loop is shown in Figure 1. For simplicity, we 
took the loop to be symmetric about its midpoint. In a ID model the loop (i.e. magnetic 
field) geometry influences the plasma evolution solely through its effect on gravity. The 
only significant restriction on the magnetic geometry is that the flux tube have a dipped 
section, as shown in Figure 1, which can support material against gravity. The geometry of 
Figure 1 would be appropriate for any of the usual prominence magnetic field models such 
as the Kippenhahn-Schluter (1957) or Kuperus-Raduu (1974), or even our recent 3D model 
( lAntiochos, Dahlburg, fc Klimchuk 1994|) . 

The loop is comprised of three distinct sections. First, is a straight vertical tube of 
height Si — this is roughly the chromospheric section. The exact position of the top of the 
chromosphere changes during the calculation. This section of the loop (s < si) is given by 


^(s) = s, hence g\\{s) = Qq, 


( 6 ) 


where z(s) is height at position s, and Qq = 2.7 x 10^ cm/s^ is the solar gravity. Next is a 


quarter circle of arc length S 2 — Si and, therefore, radius R = 2(s2 — Si)/7r. Note that the 
maximum height of the loop is Si + R. This section (si < s < S 2 ) is given by. 



(7) 


The final section is a quasi-cosine curve that has a maximum height at S 2 to match the 
quarter circle, and a minimum height at the loop midpoint, s = L, and has a total depth 






- 10 - 


D. This part (s 2 < s < L) is given by: 


/^From the form above, we note that the gravity term in the momentum equation is negative 
for s < S 2 , and positive for s > S 2 . The gravity term is continuous, but does not have 
continuous derivatives. 

The geometry parameters that we must specify (Figure 1) are the position of the top 
of the chromosphere, Si = 1 x 10® cm, the position of the loop apex S 2 = Si + 2.5 x 10® cm 
(therefore the maximum loop height si + i? = si + 5 x 10®/7r cm); the total loop half-length, 
L = 11 X 10® cm, and the depth of the coronal dip, D = 5 x lO^cm. These values are much 
more appropriate for quiescent prominence loops than the values used in our previous work, 
which are suitable only for low-lying active region prominences ( pahlburg, Antiochos, fc 
Klimchuk 19981) . 

For the radiative losses we take a simple analytic form. Since the radiative losses 
depend sensitively on elemental abundances which vary substantially on the Sun ( |(Jook ef 
|al 1989| ), there is no longer a “standard” curve to use. The important point is that our 
model does not depend on the details of the radiative loss curve; it is sensitive only to the 
spatial dependence of the heating. In the upper transition region and corona (T > 10^ 
K), we take A(T) = 1.0 x 10“^^/T. In the lower transition region (10^ > T > 30,000 K), 
A(T) = 1.0 X 10“^^T^. A temperature at the base of the model of 30,000 K is assumed. 
We would like the chromospheric region of the loop to remain at this temperature as 
much as possible; consequently, we take the radiative losses to vanish below 29,500 K. 

For 30,000 > T > 29,500 K, A(T) = 2.7 x 10-24((T - 29, 500)/500), and A(T) = 0.0 for 
T < 29, 500 K. In addition, we drop the density dependence of the radiative losses in the 
chromospheric region where the density can become extremely large due to the exponential 
increase per scale height. In the Sun’s chromosphere the radiative losses are limited by 


hence 


ttH 


s = 


-9q 


2{L — S 2 ) 


sm TT- 


S2 


S2 


)• ( 8 ) 


z(s) = Si -|- i? • 


D 

'2 


1 — cos(7r- 


S2 


S2 
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radiative transfer, which is not inclnded in onr model. As with g\\ (s) the assnmed form for 
A(T) is continnous, bnt withont continnons derivatives. 

Finally, the form of the heating needs to be specified. We set E = Eo + F^i(s, t), where 
Eq = 1.5 X 10“^ ergs/cm^/s, is the nniform backgronnd heating which stays on thronghont 
the simnlation. The spatially dependent heating Ei is ramped np over 1000 s only after 
the loop has settled into a static eqnilibrinm with the nniform heating. We set Ei to be 
constant in the chromosphere, Ei = 10“^ ergs/cm^/s, for s < Si, and have an exponential 
drop-off in the corona, 

El = 10“^exp(—(s — Si)/A), for s > si, (9) 

where the damping length A is chosen to be 10,000 km. This form for the spatial dependence 
of the heating is mnch more physical than that nsed previonsly (|Antiochos fc Klimchnk| 
1991|) , in which we assnmed that the heating was concentrated in a very narrow region 
centered at the highest point in the loop. Here, there is no special relation between the 
heating and the loop geometry. We are simply assnming that the heating propagates into 
the loop from below, and is deposited over some spatial scale small compared to the loop 
length. 


3. Results 

In order to solve Eqnations (1) - (3), appropriate bonndary and initial conditions mnst 
be specified. As bonndary conditions on the loop we assnmed a rigid, constant temperatnre 
base at s = 0, and symmetry at the top, s = L. To derive the initial eqnilibrinm we 
simply started with a discontinnons temperatnre profile in which T = 10® K in the corona 
s > Si, and T = 30, 000 K in the chromosphere. The density was constant in the corona, 
n = 0.7 X 10® cm“^ for s > si, and increased exponentially with depth in the chromosphere 
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at the appropriate scale height. Hence, the plasma was in approximate force balance to start 
with, but not in thermal equilibrium. We then let the system evolve until an equilibrium 
was reached. 

After ~ 35, 000 s the loop had settled into a static equilibrium with negligible residual 
motions, Vmax <0.1 km/s. We took this state to be the initial equilibrium, t = 0 s. 
Figures 2, 3, and 4 show the temperature, density and velocity prohles respectively at this 
time. Even though the loop has a dip geometry, the temperature and density prohles are 
exactly what one would expect for a standard coronal loop ([Rosner, Tucker, fc Vaiana] 
|1978| , [Vesecky, Antiochos, fc Underwood 1979| ). The maximum temperature occurs at the 
loop midpoint and has a value of 1.2 x 10® K. Note also that the top of the chromosphere 
has moved only slightly upward from its starting position at Si = 10, 000 km. 

We then turned on the spatially dependent heating, Ei, ramping it up linearly over 
1,000 s. The plasma responded as expected by hrst heating and evaporating chromospheric 
material into the loop before a condensation appeared at the midpoint. This evaporative 
phase of the evolution lasted a considerable time, approximately 60,000 s. Figures 2, 3, 
and 4 show the plasma prohles at three times: 10, 30, and 50 thousand seconds after the 
spatial heating turn-on. The temperature at the loop midpoint reached a maximum value 
greater than 1.7 x 10® K at around 10,000 s, but then started to decline as the radiation 
losses at the midpoint began to exceed the heating there. The position of the top of the 
chromosphere moved downward over 1,000 km as a result of the temperature increase and 
accompanying pressure increase, but then moved back to nearly its original position as the 
temperature dropped. 

We note from Figure 3 that the density throughout the loop increased steadily during 
the evaporative phase, except near the very end {t > 50, 000 s) when the condensation 
started to form at the midpoint. Due to the temperature collapse there, the pressure 
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decreased rapidly at the midpoint, which tended to evacuate the loop leg. This effect 
can be seen in the velocity profiles. At early times (~ 10, 000 s) the velocity is driven by 
evaporation, i.e., an overpressure at the top of the chromosphere, consequently, the velocity 
peaks (~ 5 km/s) near the loop base (Figure 3). The evaporative velocity decreased, 
however, as the loop temperature decreased. Later in the evolution the velocity started to 
increase again, but now peaking nearer the midpoint, which indicates that it was driven by 
the underpressure in the condensation. 

Once the condensation appeared at the loop midpoint, the subsequent evolution was 
dominated by its rapid growth. The temperature at the loop midpoint first dropped below 
the chromospheric value of 30,000 K at f = 54, 000 s. Figures 5 and 6 show the temperature 
and density profiles over the last 10,000 km of our half-loop at this time and at four later 
times in the evolution, 500, 1,000, 10,000, and 30,000 s after the condensation appearance. 
The corresponding velocity profiles over the whole loop are shown in Figure 7. 

There were two distinct phases to the condensation evolution: a transient-motion 
phase characterized by large upward and downward flows, and then a quasi-steady state 
characterized by slow upfiows throughout the loop. At first, the condensation grew very 
rapidly as the middle section of the loop dropped down to the base temperature of 30,000 
K. This temperature collapse and the resulting pressure decrease was so fast that a weak 
shock wave was launched from the loop midpoint towards the base. Figure 7. It can be 
seen from the Figure that the wave speed was 160 km/s, which agreed well with the sound 
speed in the corona. Weak shocks, like ordinary acoustic waves, are expected to travel at 
the sound speed (e.g. [ZeLdovich fc Raizer 19681) . The maximum upward velocity of the 
coronal plasma was less than 60 km/s, however, so the bulk motion was clearly subsonic 
and no strong shocks developed. The wave hit the chromosphere and then bounced back up 
the loop with a reduced amplitude. Figure 7. As shown in Fig. 7, the loop oscillated back 
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and forth for several bounces with an oscillation period ~ 1, 000 s, but after ~ 5, 000 s the 
loop settled into a quasi-steady state with a maximum upward velocity V < 5 km/s and 
with a coronal gas pressure P ~ 10“^ ergs/cm^. This flow implied a substantial enthalpy 
flux bPV/2 = 10^ ergs/cm^/s, which was sufficient to balance the excess radiative losses in 
the dipped portion of the loop, so that a steady state became possible. 

We let the simulation run for over 30,000 s after the condensation hrst appeared in 
order to verify that a true steady state was achieved. There was no signihcant change in the 
velocity prohle after the hrst ~ 5, 000 s. The evolution consisted of a very slow growth of 
the condensation. In principle, the condensation should expand until it encompasses most 
of the dipped portion of the loop in Figure 1, but this would require a time scale much 
longer than typical prominence lifetimes. The key point is that the total mass M in the 
condensation must increase as exp(z/Trg), where z is the vertical depth of the condensation 
and Hg is the gravitational scale height ~ 200 km for typical prominence temperatures. 
Hence, for a hxed velocity, the condensation’s expansion rate slows down exponentially 
rapidly with time, so we don’t expect condensations to become signihcantly larger than 
their scale heights. 


The evolution calculated by our simulation agrees very well with prominence 
observations. From the time that the midpoint temperature hrst dropped below 30,000 K, 
the condensation required 500 s to reach a half-width of 1,000 km and approximately 1,000 
s to reach a half width of 2,000 km, as shown in Figures 5 and 6. These numbers are typical 
of the time scales quoted for the formation time of prominence condensations ([Engvold| 
19761, [Zirker 1989|). Furthermore, the condensation growth slowed down abruptly after it 


reached a half-width of order 1,000 km; the half-width was only 5,000 km after 30,000 s. 
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4. Discussion 


The results described above have several interesting implications for coronal heating 
theories and for observations. The most important feature of our model is that the heating 
must be localized near the chromospheric footpoints of a loop; however, this localization is 
not severe, since we used a spatial scale of 10,000 km for the heating. On the other hand, 
our results do imply that any heating theory in which the energy input is predicted to 
occur uniformly, or localized near the loop top, would be incompatible with the observation 
of coronal condensations. Of course, even if the heating were uniform or localized near 
the top, transient condensations could still form if the heating underwent strong temporal 
variations as in the case of post-flare Hq, loops. 


Another implication of the model is that a prominence loop should be cooler and less 
dense (except in the condensation itself) than a standard, uniformly-heated coronal loop 
with the same total heating. The reason is that, in the standard coronal loop model, the 
total energy input must be balanced by the total radiation emitted {e.g., [Rosner, Tucker, 


fc Vaiana 1978|, [Vesecky, Antiochos, fc Underwood 1979|) , but in our prominence solution, 
some of the energy input goes into driving the mass flow which lifts material against gravity 
and powers the extra radiation losses at the transition region of the condensation. The 
situation is physically similar to comparing a closed coronal loop in a quiet region and an 
open fluxtube in a coronal hole. Even if the heating in each region has the same magnitude, 
the quiet-region static loop will have a substantially higher temperature and density than 
the coronal-hole steady-flow fluxtube and, hence, will appear much brighter in X-rays. 


To verify this claim we ran a simulation in which the localized heating Ei was spread 
out uniformly over the whole loop length, so that we would obtain a solution with no 
condensation. The resulting temperature and density prohles looked very similar to the 
hrst curves in Figures 2 and 3, except that the temperature maximum at the midpoint 
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had a value over 2.6 x 10® K and the density was over 2.6 x 10® cm“^. These numbers 
are to be compared to the coronal temperatures and densities of the quasi-steady state of 
the condensing loop above, 1.0 x 10® K and 1.5 x 10® cm“^ respectively. Even though the 
total heating in both cases was identical, the “coronal loop” (uniform heating case) had 
a coronal density almost two times larger than the “prominence loop” (localized heating 
case) and, therefore, would have an X-ray brightness four times larger. We conjecture 
that this effect may play a role in the origin of the dark cavities that are often seen to 
surround prominences. In our model the region of dipped magnetic flux tubes threading a 
prominence can be thought of as forming a local coronal-hole region, in which plasma is 
undergoing a quasi-steady flow into the prominence mass. Even though the low density of 
the cavity would make such flows difficult to detect, it should be informative to search for 
them with SOHO and TRACE. 

Our model predicts other flows which may be detectable by SOHO, as well. The 
key point is that if coronal heating has the type of spatial variation assumed here, then 
condensations should form in sufficiently high loops even if they don’t have dips. If the 
loop does have a dip then a quasi-steady state is possible in which the condensation settles 
to the bottom of the dip, but if the loop has no dip, we believe that no steady state is 
possible, and the loop must be constantly dynamic. We propose that this is the origin of 
the continuous dynamics observed in very high active region loops with the CDS instrument 
on SOHO ( Prekke 1997|) . Further simulations with non-dipped loop geometries are clearly 
needed. 

The main conclusion of this paper is that our model successfully reproduces the 
observed time scales for condensation formation. Not only does the condensation in the 
simulation reach a size of a few thousand km in tens of minutes, but equally important, the 
condensation stops growing signihcantly once it reaches this size. Both results are necessary 




in order to agree with observations of prominence knots. We believe, therefore, that loss 
of thermal equilibrium driven by a spatially localized coronal heating is the answer to the 
question; How does mass condense in the corona? 


This work has been supported in part by NASA and ONR. 
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Fig. 1.— Geometry of our model loop. Both horizontal and vertical distances are measured 
in units of 10,000 km. 

Fig. 2.— Temperature as a function of position along the loop for four times in the evolution. 
The temperature is measured in units of 10® K and distance in units of 10,000 km. The solid 
curve corresponds to t = 0, the time when the localized heating is turned on. The dashed 
curve corresponds to t = 10, 000 s, the dotted curve to t = 30, 000 s, and the dash-dot curve 
to t = 50, 000 s. 

Fig. 3.— As in Figure 2, but for number density in units of 10® cm“^. Note that the 
chromospheric and transition region portions of the loop are not shown so that the coronal 
density structure can be seen clearly. 

Fig. 4.— As in Figure 2, but for velocity in units of 10 km/s. 

Fig. 5.— Temperature as a function of position near the loop midpoint for hve times during 
the growth of the condensation. The temperature is measured in units of 10® K and distance 
in units of 10,000 km. The solid curve corresponds to t = to = 54, 000 s, the time that the 
condensation hrst appeared. The long-dashed curve corresponds to 500 s after to, the dotted 
curve to 1,000 s, the dash-dot curve to 10,000 s, and the short-dash curve to 30,000 s after 
to- 

Fig. 6.— As in Figure 5, but for number density in units of 10^® cm“®. 

Fig. 7.— As in Figure 5, but for velocity in units of 10 km/s and plotted over the whole 
loop. 
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